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Abstract 

We estimate the electromagnetic effect of the spectator charge on 
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\ the momentum spectra of tt + and ir produced in peripheral Pb+Pb 

collisions at SPS energies. We find that the effect is large and results 
in strongly varying structures in the xj? dependence of the ir + /ir~ 
! ratio, especially at low transverse momenta where a deep valley in the 

above ratio is predicted at xf ~ 0.15 - 0.20. 

It appears that the effect depends on initial conditions. Thus, 
it provides new information on the space and time evolution of the 
CJ» ' non-perturbative pion creation process. 
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PACS: 25.75.-q, 12.38.Mh 



^ 1 Motivation 



The process of pion production in high energy hadronic reactions (like p+p, 
p+nucleus or nucleus+nucleus collisions) belongs to the domain of the strong 
interaction theory: quantum chromodynamics (QCD). For an enormous ma- 
jority of pions produced, the mechanism underlying this process is non- 
perturbative. Here, QCD does not lead to quantitative predictions. Instead, 
phenomenological models are used [[]. These contain arbitrary scenarios of 
the hadronic collision dynamics, and evidently rely on experimental input to 
differentiate right scenarios from wrong ones. 

Most of the experimental input collected so far (at ISR, SPS, RHIC, 
and others) provides exclusively information about final state momentum 
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space (p x , p y , p z ). Observables such as inclusive distributions of kinematical 
variables, or their correlations, give no information on how the reaction occurs 
in position space (x,y,z) or in time (t). To get such information, a specific 
phenomenon directly dependent on the reaction evolution in space and time 
has to be isolated. The HBT effect 2J can be quoted as an example of such 
a phenomenon. 

This paper considers another possible candidate for such a phenomenon 
in ultrarelativistic heavy-ion collisions. We mean here the electromagnetic 
interaction between the remnants of the two nuclei (spectator systems) and 
the positive and negative pions produced in the course of the collision. The 
two highly charged spectator systems moving at relativistic velocities gener- 
ate a rapidly changing electromagnetic field which modifies the original pion 
trajectories. This causes a distorsion of observed kinematical pion spectra. 
It is to be expected that this distorsion is interrelated to the dynamics of the 
collision, and in particular to the time evolution and initial conditions of the 
participant and spectator zones. 

Our aim is to study this electromagnetic "Coulomb" effect for the specific 
case of peripheral Pb+Pb collisions at SPS energies (158 GeV/nucleon beam 
energy, ^^^=17 GeV). We choose peripheral collisions because they are 
characterised by the largest spectator charge. Our specific tasks are: 

1. to provide a general description of main features of the effect in the 
framework of a simplified but realistic, relativistic two-spectator model. 

2. to define the regions of produced pion phase-space that are particularly 
sensitive to the electromagnetic influence of the two spectator systems. 

3. to get a first idea on whether the electromagnetic effect depends on the 
evolution of the pion production process in space and time. 

Seen in a more general context, our analysis is complementary to various 
studies of the role of specific phenomena like multiple collisions 3J, isospin 
effects the neutron halo 5j, or Fermi motion |B] in particle production in 
heavy-ion reactions. To understand the role played by all these effects seems 
to us a necessity if a realisitic description of the non-perturbative heavy-ion 
collision dynamics is to be provided. 

The remainder of this paper is organized as follows. The general context 
of Coulomb interactions in nuclear collisions is shortly discussed in Section 2. 
The general idea and the details of our model are described in Section 3. The 
results of our study are presented in Section 4. Conclusions are formulated 
in Section 5. 
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2 Coulomb Effects at Low and High Energies 



Various types of Coulomb effects were investigated both experimentally and 
theoretically in the past; some examples will be enumerated below. 

Various results are available in the low (or intermediate) energy regime, 
up to a few GeV/nucleon. On the experimental side, it was realized in studies 
at Bevalac that the Coulumbic field causes extremely large effects at forward 
angles, i.e. at projectile velocity 0IBj- The Coulombic effect was also ob- 
served at subthresold energies [Hj. In studies of p+Au and He+Au reactions 
at several GeV/nucleon, the non-relativistic approach to the Coulomb field 
brought precise information on the evolution of the nuclear fragmentation 
process jTU]. Finally, an interesting charge asymmetry was found in the ki- 
netic energy spectrum of particles produced in collisions of 9 GeV protons 
with emulsion nuclei jllj . 

The situation on the theoretical side is rather complicated. Attempts 
were made to describe the intermediate energy minimum bias data by a 
multiparametric modelling of the pion source fireballs, and calculating the 
subsequent transport of pions in the Coulomb field of the spectators [TJl HHj . 
Quantal effects were discussed in [Hj. Other simple models were considered 

in nana Ed. 

In the high energy regime, studies of spectator Coulomb interactions offer 
two advantages relative to the low energies. On the experimental side, im- 
pact parameter (centrality) selection becomes possible by produced particle 
multiplicity measurement or by forward calorimetry. This allows a selection 
of peripheral collisions where the spectator effect will be the largest. On the 
theoretical side, seen in the context of typical relativistic pion velocities, the 
break-up of the spectator system can be assumed to be slow. This simplifies 
the picture. However, a fully relativistic description of electromagnetic effects 
involving e.g. retardation phenomena becomes a necessity at high energies. 

In early cosmic ray experiments, large tt + /tt~ asymmetries were observed 
for produced slow pions [TH]. As far as nucleus+nucleus collisions and SPS 
energies are concerned, various experimental fHJ E] arid theoretical 
papers on the influence of the electromagnetic field on particle production in 
the mid-rapidity (xp — 0) region 1 of Pb+Pb reactions are known. However, 
these focus on the role played by the electromagnetic field originating from 
particles taking direct part in the reaction. 

As far as the influence of the two spectator systems in the rest of avail- 
able phase-space is concerned, we are aware of only one published measure- 

1 Note: the Feynman variable xf = ^Pl/^/s, the transverse momentum px, the energy 
E, and other related quantitites will always be considered in the nucleon+nucleon c.m.s. 
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ment j2Hj- This has been performed for a very limited acceptance range 
(forward angles i.e. Pt — 0). New, better measurements can be expected in 
future from the NA49 apparatus j23]. No theoretical description is known to 
us. To provide such a simplified but realistic description is the main purpose 
of the present paper. 



3 Modelling a Peripheral Pb+Pb Collision at 
SPS Energies 



Figure 1: Our simplified view of a Pb+Pb collision. The hypothetical pion 
emission region is reduced to a single point in position space (see text). 

Our aim is to obtain a realistic estimate of the spectator-induced elec- 
tromagnetic effect, avoiding however the discussion of complex and poorly 
known details of soft hadronic particle production. Thus we decide on a 
maximally simplified approach, illustrated in Fig. ^ A peripheral Pb+Pb 
reaction can be imagined to consist of three steps: 

1. the collision takes place at a given impact parameter h. The two highly 
charged spectator systems follow their initial path with essentially un- 
changed momenta; 

2. the participating system evolves until it finally gives birth to final state 
pions. The evolution of the pion emission region in space and time is 
a priori unknown; 

3. charged pion trajectories are modified by electromagnetic interaction 
with the spectator charge; the spectator systems undergo a complicated 
nuclear fragmentation process. 




^ projectile 
+ + \ spectator 



b 



target ^ 
spectator 
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We model these steps in the following simplified way: 

1. a peripheral Pb+Pb collision involving 60 participating nucleons is as- 
sumed. This corresponds to an impact parameter b equal to 10.61 fm 
(see Sec. 13. lj) . The two spectator systems are modelled as two uniform 
spheres in their respective rest frames. The sphere density is the stan- 
dard nuclear density p = 0.17/fm 3 . The total positive charge of each 
spectator system is Q = 70 elementary units. In the collision rest frame 
the two spheres become charged disks (Fig. [TJ). 

2. the pion emission region is reduced to a single point in space, namely 
the original interaction point. The emission time £e is a free parameter 
in our model. For peripheral Pb+Pb reactions studied here, we assume 
that the initial two-dimensional (xf,Pt) distribution of the emitted 
pion is similar to that in underlying nucleon+nucleon collisions. 

3. charged pions, with their initial momentum vector defined above in 
point 2., are numerically traced in the electromagnetic field of the spec- 
tator charges until they reach a distance of 10,000 fm away from the 
original interaction point and from each of the two spectator systems. 
The fragmentation of the spectator systems is neglected; the influence 
of participant charge, strong final state interactions, etc, are not con- 
sidered. 

For each of the above steps, the details of our approach are explained 
below. 

3.1 Collision Geometry 

We adjust the geometry (centrality) of the considered Pb+Pb collision to 60 
participating nucleons in order to make it comparable to that of peripheral 
heavy-ion data samples collected at the SPS j2Sj- 

The relation between the reaction impact parameter b, the number of 
participating nucleons N part and the spectator charge Q is defined by the 
nuclear density profile and the elementary nucleon+nucleon cross section. We 
study this relation by means of a geometrical Monte Carlo simulation. This 
simulation (Fig. EJ) produces spatial distributions of protons and neutrons 
using nuclear density profiles obtained for 208 Pb by Mizutori et al. in the 
Hartree-Fock-Bogoliubov (HFB) approach For each Pb+Pb reaction, 

a nucleon is defined as participant if it is crossed by one or more nucleons 
from the adverse nucleus within a transverse radius of less than 1 fm. This 
corresponds to an elementary nucleon+nucleon cross-section of 31.4 mb in 
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Figure 2: Geometrical model of a Pb+Pb collision. The filled ellipses repre- 
sent participating nucleons while empty ellipses represent spectator nucleons. 



good agreement with experimental data [21] ■ No re-interaction with the 
spectator system is considered. 

Technically, the simulation is performed at a fixed value of the impact 
parameter b. We obtain the average number of participants N part = 60 
at b = 10.61 fm. At the same time we obtain 2 the average charge of the 
spectator system Q « 70 e.u. 

Another important geometry parameter to be estimated is the displace- 
ment Ab of the spectator protons' center of gravity relative to the center of 
gravity of the original Pb nucleus. This displacement originates from trivial 
geometrical reasons (Fig. |3J). Our simulation gives Ab = 0.76 fm. Thus the 
effective distance of closest approach between the centers of gravity of the 
two spectator systems will be b = b + 2Ab = 12.13 fm. 

We consider that the exact shape of the spectator systems is not essential 
for our subsequent simulations of electromagnetic effects. For the sake of 
clarity, we model the two spectator systems as two homogenous spheres with 
standard nuclear densities p=0.17/fm 3 and with properly shifted centers of 
gravity (Fig. [TJ). As discussed above, the distance of closest approach of the 
two spheres' centers of gravity is b = 12.13 fm. The total charge of each 
sphere is Q = 70 elementary units. 

For the case of peripheral collisions considered here, estimates on geom- 
etry parameters may strongly rely on a proper description of the nuclear 
profile. Therefore, in order to cross-check our results, we performed an in- 
dependent simulation using a different nuclear profile: a phenomenological 
parametrization obtained by Trzcihska et al. on the basis of antiprotonic 

2 The spectator charge Q is a numerical result of the simulation. It does not a priori 
follow from N part due to possible neutron halo effects |2Hj ■ Our Monte-Carlo takes account 
of these effects. 
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Figure 3: As some of the nucleons become participants, the center of gravity 
of the spectator protons is displaced relative to the center of gravity of the 
original Pb nucleus. 

atom data With respect to the Mizutori profile, the deviations on b' , Q 
and Npart are small and remain below 1.5%. Thus we consider our estimates 
as fairly precise. 

3.2 Initial Pion Emission 

Having established the initial geometry of the peripheral Pb+Pb reaction, 
we subsequently model the emission of produced 7r + and n~ . 

We reduce the unknown initial emission region to a unique point in space, 
namely the original interaction point (Fig. HJ. We assume one emission time 
t E which is a free parameter of our model; various values of t E will be con- 
sidered. Such a simplification of initial conditions gives a convenient way 
to estimate the dependence of the electromagnetic effect on the characteris- 
tics of the initial pion emission process (like the pion formation time or its 
distance from the two spectator systems). 

We assume that the initial kinematical spectra of the emitted pions are 
similar to these in nucleon+nucleon collisions and that they follow wounded 
nucleon scaling [SU]. Full azimuthal symmetry of the emission is assumed. 
The detailed shape of the emitted pion density distribution requires an ad- 
ditional discussion. 

A very precise and complete set of experimental data on ir + and 7r~ pro- 
duction in p+p collisions is available at the SPS [27|. The published double 
differential spectra of pions in xp and px display a strong dependence on pion 
isospin. They also contain local shape structures attributed to production 
of hadronic resonances and their subsequent decay into pions. It has been 
demonstrated that such elementary effects reappear in nuclear collisions [3]. 




Figure 4: Assumed invariant density of tt + emitted per participant pair in 
Pb+Pb collisions, drawn as function of xp at fixed pt values. From top to 
bottom, the curves correspond to px = 50, 100, 200, 400, 600, 800, and 
1000 MeV/c. For clarity, the two top curves are multiplied by 2.5 and 1.5, 
respectively. The invariant density of emitted tt~ is assumed to be identical 
to the above (see text). 



Any attempt at an in-depth description of Pb+Pb interactions must take 
them into account. 

Such a description is, however, beyond our specific scope. As we con- 
centrate on the spectator-induced Coulomb effect per se, we need a clear 
distinction between hadronic and electromagnetic phenomena to allow for 
an easy interpretation of our results. Thus, we simplify the situation: 

1. we neglect isospin effects i.e. assume equal initial emission spectra for 
7r + and 7r~; 

2. we construct a smooth two-dimensionnal shape that reproduces the 
most basic features of the pion production data j2Zj; it does not include 
the more subtle, local shape structures. 
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The resulting two-dimensional invariant pion density 3 is drawn in Fig. 
We do not consider values higher than \xp\ = 0.8. We checked that below 
Pt = 1 GeV/c, our curves typically deviate from measured average pion 
( - y~ ) densities by about 10%. Maximal deviations reach up to 21% of 
our assumed density. Such a gross description is sufficient for our specific 
aim. Evidently it cannot be used as a "fit" nor "parametrization" of the 
data. Note that the failure of oversimplified analytical parametrizations to 
reproduce the p+p data has been demonstrated in [22] . 

For completeness, we write the analytical form of our emitted pion density 
per participant pair: 



2 d?N 

= ^ a n exp(-(x/b n ) Cn )exp(-u T /d n ) , (1) 



dp" pb+Pb-,TrX „ 12 



where it = ti + or tc , N part = 60 (Sec. 13. 1J) . x = \J x 2 F + g 2 , ut = \J q 2 + Pt, 
and numerical parameter values are listed below: 



n 


a n [c 3 /GeV 2 ] 






d n [GeV/c] 




1 


2.32229 


0.369967 


2. 


0.191506 


g= 0.01 


2 


24.4563 


0.0873833 


1.001 


0.12 


q = 0.334968 [GeV/c] 



3.3 Propagation of Charged Pions in the Spectator 
Electromagnetic Field 

The initially produced charged pions are subjected to the electromagnetic 
field of the two spectator systems moving at relativistic velocities. 

We choose the overall center of mass system of the collision to calculate 
the evolution of pion trajectories. For symmetric Pb+Pb collisions and ne- 
glecting the Fermi motion effect (see this is also the nucleon+nucleon 
center of mass system. 

We assume that for ultrarelativistic Pb+Pb reactions, the nuclear frag- 
mentation process is slow relative to the relativistic pion velocities and can 
be neglected in a first approximation. Thus the spectator nucleon velocity 
remains constant and identical to the velocity of the parent Pb ion. 

We define our space coordinate system in the way shown in Fig. For 
simplicity, we label the two spectator systems as left (L) and right (R). 
We set the time scale such that at t — the center of gravity of each of 

3 Note that in terms of xf, Pt and azimuthal angle cj), our invariant pion density per 
event and per participant pair writes ^ = ^ p ^E_ _^f|L_ . 



9 




(0,0,0) 1 



target 
spectator (L) 



projectile 
spectator (R) 



n 



Figure 5: The space coordinate system for the Pb+Pb collision. The origin 
(x, y, z) = (0, 0, 0) is located at the interaction point. The projectile and 
target spectators are denoted left (L) and right (R). We define the impact 
parameter b as a vector. 



the spectator systems is found at zl = zr = 0. Thus the time-dependent 
position of the moving spectators is: 



R L {t) = -b/2 + v L -t, 

R R {t) = b/2 + v R -t. (2) 

Here, Vl and Vr are the two spectator velocity vectors. For the SPS energy 
discussed here vl = vr = vs = 0.994c. The value of the impact parameter 
vector b is set to b = 12.13 fm (see Sec. 13. 1|) . 

Each of the two spectator systems defines its own rest frame. We define 
E' L as the constant electrostatic field generated by the spectator L in its rest 
frame, and E" R as the field generated by the spectator R in its rest frame. 
Having assumed both systems as uniform spheres with a normal nuclear 
density p = 0.17/fm 3 and with a total charge Q = 70 elementary units 
(Sec. 13.1 J) we write: 

E> (p) = { k ® rf for r' c > R s 
^ L ' I k Q r'J R? for r' c <R s 1 ' 

E"(r») = l kQ r "' /r " 3 fOT r " >Rs U) 
' ' k Q P>/R s 3 for r" c <R S K J 

In the equations above, r' c (r") is the position relative to the center of the L 
(R) spectator, defined in this spectator's rest frame, k ~ 1.44 MeV • fm/e 2 
is the electrostatic constant. R s = [A^ spec /(4/37rp)] 1//3 is the sphere radius 
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defined by the number of spectator nucleons N spec . For the case considered 
here R s = 6.3 fm. 

We transform the fields E' L , E R to the center of mass system. Here, the 
moving spectator charge generates both electric and magnetic fields. From 
the general Lorentz transformation [HI] we get 



E L (r, t) = ls E' L (rl) - ^ ( ^ • E' L {r' c , 

7 S + 1 c V c 

B L if,i) = 7 S x E' L {r'^j 



cl I 1 



(5) 



for the left spectator and 



E R (r, t) = 7s £» (r» ) - ?* ( ^ • P R {? ) 



(6) 



for the right spectator. 

In the equations above, the j s factor is defined as 7s = (1 — t^/c 2 ) -1 / 2 . 
The vectors E L (E R ) and B L (B R ) are respectively the electric and magnetic 
fields generated by the left (right) spectator at the space-time position (r, t). 
The relation between r = (x, y, z) and the spectator rest frame coordinates 
r' c = (x' c , y' c , z' c ) , r" = {x" c ,y" c -,z" c ) is given by the Lorentz transformation and 
the impact parameter vector b = (b x ,b y ,0): 

x' c = x + b x /2 ( x" c = x - b x /2 

y' c = y + b y /2 | y'> = y-b y /2 (7) 

z' c = ^ s {z + v s t) ( z" c = 7 s (z - v s t) 

We now consider a charged pion emitted at time t = tg from the inter- 
action point r = (0,0,0) with its initial momentum p n (t = tg) specified by 
Eq. from Sec. 13.21 This situation is illustrated in Fig. El The Lorentz 
force acting on the pion is: 

^ = F„(r, t) = q n fE(r, t) + x B(r, t) 

where q n is the pion charge, E(f,t) = E L (f,t) + E R {r,t) and B(r,t) = 
Bi(r, t) + B R (r, t) are standard superpositions of fields from the two sources. 
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The resulting pion trajectory f n (t) is denned by its time-dependent ve- 
locity v n (r, t): 

d 4 = ^t) = -^£=, (9) 

where m n is the pion mass. 

An important feature of Eq. (jHJ) is that it implicitly takes account of rel- 
ativistic retardation effects. This feature has been explicitly demonstrated 
in |3I]. 

Technically, the propagation of the pion is made by means of an iterative 
Monte-Carlo procedure. This procedure starts at f — (0,0,0) and t = t E 
and calculates the Lorentz force F n (f,t) and the corresponding change of 
pion momentum and position in small steps in time. The variable step size 
depends on the actual distance of the pion from the nearest spectator system. 
The procedure is iterated numerically until the distance of the pion from 
the origin (0, 0, 0) is r > R ma x and at the same time, the distances of the 
pion from the spectators in their respective rest frames are simultanously 
r' c > R max and r" > R ma x- We found that for geometries considered in the 
present paper, the value of the parameter R max = 10,000 fm is sufficiently 
large to reproduce asymptotic momenta. The procedure is weighted, that is, 
each pion is generated with its proper weig ht deduced from Eq. (0). 

After the propagation procedure is finished, the same weight is used to fill 
the final state pion spectra. Negatively charged pions that do not escape 
from the spectator potential well are rejected by our procedure and do not 
enter into the final state distribution. 



4 Results 

In this Section the results of our Monte-Carlo studies are discussed. 



4.1 Charged Pion Spectra 

We start by estimating the influence which the spectator charge exerts on 
double-differential spectra of pions produced in peripheral Pb+Pb reactions. 
For the sake of clarity, this part of the discussion remains limited to the 
simplest situation where the pion emission time t E is equal to zero (immediate 
pion creation). This description will be generalized in Sec. I4.3l bv considering 
different t E values. 

The situation is illustrated in Fig. El Panel (a) shows the initial spectra of 
emitted pions. As explained in Sec. I<1 2| in our simple model these spectra are 
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Figure 6: Double-differential density of positively and negatively charged pi- 
ons, produced per participant pair in peripheral Pb+Pb reactions, a) Initial 
density of emitted 7r + and 7r _ . b) Density of tt + in the final state of the 
reaction, c) Density of ir~ in the final state of the reaction. In all the panels, 
the pion density is drawn as a function of xp at pp = 25 MeV/c (thin solid), 
75 MeV/c (dash), 125 MeV/c (dot), 175 MeV/c (dash-dot), and 325 MeV/c 
(thick solid); each of these pr values corresponds to a bin of ±25 MeV/c. 
This simulation was made assuming the pion emission time equal to zero. 
At negative xp reflected curves are drawn. 
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identical for n + and n . The presented , d ^ density distributions (scaled 
down by the number of participant pairs) follow directly from Eq. (0). 

In panel (b), the corresponding distributions of tt + in the final state of 
the Pb+Pb reaction are shown. These are obtained by our Monte-Carlo 
simulation described in Sec. 13.31 It is clearly apparent that the distribu- 
tions are distorted by the Coulomb repulsion between the pion and spectator 
charges. The effect is largest for pions moving close to spectator velocities 
(xp ~ ±0.15) and at low transverse momenta (pp = 25 MeV/c). Here, two 
deep valleys in the 7r + density are visible. A similar but smaller distorsion is 
also apparent at pp = 75 MeV/c. 

An opposite distorsion is present for 7r~ densities shown in panel (c). 
Negative pions are attracted by the positive spectator charge and gather at 
low transverse momenta close to spectator velocities. This results in the 
presence of two large peaks at xp ~ ±0.15. Remnants of these peaks are 
apparent at pp = 75 MeV/c. 

Thus, the Coulomb field induced by the spectator charge at SPS energies 
appears strong enough to produce visible distorsions on the pion d fJ^ pT dis- 
tribution, which extends over several orders of magnitude in terms of pion 
density and over many GeV in terms of kinetic energy 4 . 

4.2 7r + /7r~ Ratios 

Further, more precise information on the spectator electromagnetic effect 
can be obtained by considering the density ratios of produced positive over 
negative pions: 7i + /tt^. In our model, the 7r + /tt" ratio for initially emitted 
pions is by definition equal to unity in the whole phase space (Sec. 13. 2j) . Any 
deviation of this ratio from unity in the final state of the Pb±Pb reaction is 
therefore directly resulting from the spectator-induced Coulomb interaction. 

Fig. [3(a) shows the x^-dependence of the n + /n~ ratios computed for the 
case considered in the preceding Section, namely the pion emission time tp 
assumed equal to zero. The spectator Coulomb field appears now to produce 
a characteristic, complex pattern of deviations from unity. The first element 
of this pattern is a double, two-dimensional valley which covers the low- 
Pt region in the vicinity of xp ~ ±0.15; the valley remains still visible at 
Pt = 175 MeV/c. The second element is a smooth rise of the 7t + /tt~ ratio 
at higher \xp\. This rise is present for all the considered pp values; at fixed 
xp, the ratio slowly decreases with increasing p T . 

4 In this context, it is interesting to note that the initial spectator c.m.s. electrostatic 
potentials induced on the pion are relatively small (-jp ~ 17 MeV for the case considered 
in Fig. EJ). 
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Figure 7: Ratio of density of produced 7r + over produced 7r~ in the final 
state of the peripheral Pb+Pb reaction, obtained for five values of the pion 
emission time tp- In all the panels, the 71 + /tt~ ratio is drawn as a function 
of x F at p T = 25 MeV/c (thin solid), 75 MeV/c (dash), 125 MeV/c (dot), 
175 MeV/c (dash-dot), and 325 MeV/c (thick solid). The small cusps on 
the curves correspond to the statistical fluctations in our Monte Carlo. At 
negative x F reflected curves are drawn. 
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Figure 8: Final state 7i + /tt~ ratio for the peripheral Pb+Pb reaction, shown 
at fixed values of pr as a function of xf- The five considered values of the 
pion emission time t# are differentiated by means of different line types. At 
negative xp reflected curves are drawn. 
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4.3 Dependence on Initial Conditions 

The central issue of this paper is the sensitivity of spectator-induced Coulomb 
effects to initial conditions of the pion production process. This issue is stud- 
ied in Fig. El(a-e) where the 7r + /7r~ pattern obtained with the pion emission 
time t E equal to zero is compared to these computed with t E = 0.5, 1, 1.5, 
and 2 fm/c. 

Indeed, clear differences appear for different t E values. With increasing 
t E , a broadening of the double valley at x E ~ ±0.15, and a decrease of the 
7r + /7r~ ratio at higher absolute xf, are visible. 

The study is further quantified in Fig. |HJ Here, the n + /7T~ patterns ob- 
tained with different t E values are directly compared at fixed p?- It becomes 
evident that a change of 0.5 fm/c in the pion emission time is sufficient to 
produce visible changes in the it + /tt~ ratio. In the low transverse momentum 
region (pt < 200 MeV/c), a displacement of the double valley towards higher 
absolute x E becomes apparent on top of its broadening with t E . At higher 
transverse momenta, the ratio seems to stabilize for t E > 1 fm/c. 

This evident sensitivity of the 7t + /tt~ ratio to the pion emission time 
t E clearly indicates that the Coulomb effect induced by spectator protons 
depends on initial conditions imposed on pion production. Let us remind that 
the different considered values of t E are, in our model, equivalent to different 
distances between the pion formation zone and the two spectator systems. 
It can therefore be concluded that the Coulomb distorsion pattern induced 
by spectators on ir + and 7T~ spectra carries information on the evolution of 
the particle production process both in space and in time. 

5 Summary and Conclusions 

The electromagnetic interaction between the spectator protons and charged 
pions produced in the peripheral Pb+Pb reaction has been studied by means 
of a simplified model. This interaction is strong enough to produce visible 
distorsions in the final state densities of positive and negative pions. 

The main feature of this "Coulomb" effect is a big dip in the 7r + density 
distribution at low transverse momenta in the vicinity of x E ~ ±0.15, ac- 
companied by an increase of 7r~ density in the corresponding region of phase 
space. This results in the presence of a double, two-dimensionnal valley in 
the 7r + /7r~ density ratio. At higher absolute x E another distorsion, namely 
a smooth increase of tt + /tt^ with x E , may appear. 

The sensitivity of this electromagnetic effect to initial conditions imposed 
on pion production has been estimated. The effect is clearly sensitive to 
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initial conditions. Changes of the pion emission time by 0.5 fm/c (in c.m.s. 
time) are sufficient to modify the observed distorsion pattern. In our model, 
such changes are equivalent to changes of position of the formation zone by 
0.5 fm relative to the two spectator systems. Thus, the electromagnetic effect 
appears to depend on the evolution of the pion production process in space 
and in time. 

Our study demonstrates the importance of new, double-differential data 
on the xf and pr-dependence of pion production in peripheral nucleus+ nucleus 
collisions. Further analyses of this subject should take account of other effects 
present in nuclear reactions (multiple collisions, isospin, nuclear fragmenta- 
tion and possibly strong final state interactions of pions with the spectator 
systems). It seems nevertheless clear that the electromagnetic phenomena 
induced by the presence of spectator charge carry interesting information on 
the mechanism of the non-perturbative particle production process. 
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